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ABSTRACT 

In the framework of a field research of the dynamics and renewal of water in Montevideo’s Bay, 3.7 x10? Bq of 
tritiated water were distributed in the north-east region of the bay. Remnant total tritium activity in the bay was 
estimated from the activities measured in 20 volume cells during 100 hours. The time course of the tail of the 
tritium remnant function suggests an asymptotic non-exponential water renewal process. Here a previous 
conventional two compartments model is generalized by the introduction of a third fractal compartment. A 
mathematical model compatible with both, fractional and fractal compartments analysis, is constructed. A 
reduced order model is derived and a general analytical solution to describe tritium’s kinetics is obtained by 
perturbation methods. A detailed solution in the framework of fractal compartment analysis is constructed, a 
formula to estimate renewal time in the bay is derived, and an asymptotic formula for tritium remnant function is 
given. The introduction of the fractal compartment allows an explanation to the asymptotic behavior of measured 
tritium’s remnant function. The obtained formulae can be applied to the non-exponential wash out kinetics from 
the bay related with the trap effect characteristic of the fractal compartment. This could be of interest to follow 
and assess the impact of the liberation in coastal waters of toxic chemical substances or dangerous microbial 


populations, after a sudden contamination accident in the harbor. 


Key words: fractal compartments, fractional compartments, mathematical modeling, radioactive tracers, surface 


hydrology, tritium. 


INTRODUCTION 

This draft paper describes the construction of a new three compartments model for water renewal in 
Montevideo’s Bay that improves a conventional two compartments model introducing a fractal 
compartment. 

Montevideo’s Bay is located in the left (north) bank of River Plate estuary, in the transition region 


where the salt wedge ends (Figure 1). 
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Figure 1: A view of Montevideo’s Bay and adjacent coastal waters. 


Semidiurnal astronomical tides, winds and river discharges produce variable effects both on coastal 
water dynamics adjacent to the bay and on water dynamics inside the bay. 

The bay has an average surface area of 1.6 x10’ m? and an average volume of 4.8 x10’ m°, with an 
average depth of 3 m. Although water movements in this fairly shallow water body are highly 
influenced by the wind, certain regular patterns emerge related both with tidal motions and permanent 
vertical stratification of bay’s waters. 

During rising tides the average movement of coastal waters is toward the West and usually fills the bay 
by a mechanism of displacement of contiguous water masses. During decreasing tides the water 
moves towards the East, produces a clear-cut circulatory flow in clockwise direction and the bay is 
partially emptied. 

According to the available evidence, the highly polluted Montevideo’s bay operates as an intermittent 
tidal pump that injects contaminants in the nearby coastal waters, with potential effects on the beaches 
located east to the bay’s mouth [1]. As consequence, it could be of certain importance to quantify the 
wash out kinetics of toxic chemical substances or dangerous microbial populations, introduced in the 


bay after a sudden contamination accident in the harbour. 
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In order to assess this effect, several field studies were done in the bay itself and in the adjacent coastal 
waters covering water and sediments dynamics [2]. In one of these studies, 3.7 x10’? Bq (100 Ci) of 
tritiated water (THO) were evenly distributed in the north-east region of the bay, during five hours of 
continuous injection from a moving boat. Figure 2 shows a sketch of the bay and adjacent coastal 


waters. The region where the tritium was uniformly dispersed is shown in light gray. 
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Figure 2: Montevideo’s bay and proximal coastal waters. The main injection points of polluted fresh 
waters and suspended solids in the bay in the time when the tritium experiment was done are indicated by 


the white arrows. 


The whole bay was divided in 20 concentration cells, applying a methodology analogous to the one 
used to elaborate rainfall charts, taking into account available bathymetric data and corrections from 
field data obtained in situ during field work. Tritium concentrations (activities per unit volume) and 
other parameters (temperature, electrical conductivity, etc.) were measured in vertical profiles during 
three weeks, in the mid-point of each cell, first twice a day and the on a daily basis. Remnant total 


tritium activity was estimated from cells volumes and midpoint cells activity concentrations. 
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A best fit straight line was determined from experimental data, excluding the first 18 hours since the 
injection. The slope allowed estimate a global renewal time of 29 hours. This result was interpreted (in 
a coarse approximation) in the framework of a one compartment model. However the extrapolation to 
the time of injection gives a tracer value higher than 100 %. This fact jointly with the spatial and 
temporal details of the measured tritium kinetics, available bathymetric data, information about water 
movements in a tidal environment measured with drogues, fluorescent tracers and current meters, as 
well as the results of computer fluid dynamics simulations (using numerical models averaged in 
depth), suggests that (in a less coarse approximation) the bay can be meaningfully divided in two main 
compartments [3]: a North-East (NE, average volume 1,6x10’ m*) and a South-West (SW, average 
volume 3,2x10’ m?) compartment. So a linear classical two compartment model was constructed and 
analysed, with tritium activities in the SW and the NE compartments as state variables [4]. From now 
on we will call it the old model. 

However, the time course of the tail of the measured tritium total remnant function and of the tails of 
the measured activities in the NE and SW compartments (long tails of activity that seems to vanish 
slower than an exponential) points to an effect of anomalous dispersion of THO, perhaps related with 
constrained random walks in some places of the system. 

The distribution of dead zones jointly with details of the geometry of the bay bottom, including the 
distribution of a porous medium formed by fluid mud, suggest the introduction of a third fractal 
compartment, besides the original classical compartments (SW and NE). 

In the framework of nonconventional compartments kinetics there are two approaches to describe the 
exchange tritium flows between a fractal compartment and the classical ones. One approach is known 
as fractal compartment analysis [5]. In linear fractal compartment analysis the kinetic coefficients that 
are multiplied by the concentrations to obtain the flows to or from a fractal compartment, are 
represented as decreasing functions of time that vanish as power laws. The other approach is called 
fractional compartments analysis [6]. In linear fractional compartments analysis the flows to or from a 
fractal compartment are proportional to fractional derivatives of the concentrations. Both approaches 
allow us to take into account constrained transport processes in a heterogeneous medium like the 


fractal compartment, but they are not fully equivalent. 
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The purpose of the present paper is threefold: 

(1) To generalise the old model to a linear three compartment model, in the framework of an extension 
of compartments theory that allows taking into account anomalous dispersion in mass exchange 
processes involving a fractal compartment, and is compatible with both, fractal and fractional 
compartments theories. 

(2) To apply the asymptotic methods of perturbation theory to a reduced order model derived from the 
linear three compartment model, obtaining an approximate general analytical solution for the activities 
in each compartment and for the renewal time, compatible with both fractal and fractional 
compartment analysis. 


(3) To derive an asymptotic formula for the time course of the tail of the tritium remnant function. 


A GENERALISED THREE COMPARTMENTS MODEL FOR TRITIUM KINETICS 

To take into account possible anomalous diffusion effects, we propose a model based on three 
compartments: the two NE and SW of the old model and a third one F, fractal. This seems to be the 
simplest model structure that allows a description of the slow non-exponential vanishing of the trittum 
activity suggested by the experimental results, while retaining the possibility to describe the initial 
tritium kinetics already achieved by the old two compartments model. Figure 3 shows a sketch of the 


proposed three compartments model. 


Figure 3: Sketch of the three compartment model. 
In Figure 3 the water flows exchanged between the NE and SW compartments and between the coastal 
waters and the SW compartment, have two subscripts. Here we use the engineering convention to label 


5 
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the exchange volume or mass flows between compartments: the second (right) subscript indicates the 
compartment from which the flow comes from and the first (left) subscript indicates the compartment 
to which the flow is directed to. 


Volume balance: 


Let us consider the bay split in a SW compartment of volume V,(t) , a NE compartment of volume 
Vi (t ) and a fractal compartment (F) of volume V, (t ) so that the bay’s volume is 


V(t) = V(t )+V,(t )+V, (t ) (where f represents the considered time instant, as usual). From now on 


we suppose that volume flows between the fractal compartment and the other two compartments are 
negligible compared with the volume flows exchanged between the NE compartment, the SW 


compartment and coastal waters. 


The SW compartment (compartment 1) receives three volume inflows: one inflow Q,, (t ) from 
adjacent coastal waters, through the bay’s mouth, another inflow Q, at ) from the NE compartment 


and a third inflow QO alt ) of sewage waters injected from the land into the harbour’s docks (these 


docks are located in the SW compartment). This same compartment sends two volume outflows, one 


0; (t ) to the adjacent coastal waters and the other Q,, (t ) to the North East compartment. The 


volume balance for the SW compartment is given by: 


ZV )= Qc} Q.ilt)-Onlt)+ Galt) Onl) 


The NE compartment (compartment 2) receives a fresh volume inflow O,,(t ) (due to Pantanoso 


creek, Miguelete creek, as well as the 12 de Diciembre and other sewage injections from land shown 


in Figure 2). Its volume balance is given by: 


(2) <Vj(t)=Qnlt)-Qrlt)+ Opal) 


Activity balance: 
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Ifa, (t)= Ait) , a,(t)= Axle) and a; (r) = an represent the tritium activity concentrations 
3 


interpretable as averaged values in the SW, the NE and the F compartments, respectively, neglecting 
the backflow of activity from the adjacent coastal waters, the activity balances for each compartment 


in the proposed model are: 


(3a) # 4 ()=0.4()+0,,(0))- tt + Qra(t)- a és) 219) | il i 


(3b) —A,(t)=@,,(t)- Ate 9.(0 


Go) Lae) 8 A) anf |e rae 


V(t) V3(t) V5 (t) 


2 | Aalt) 

—~K,,| 32 

eg 
While the NE and SW compartments are considered as the homogeneous kinetic spaces of classical 


compartment theory [5], the third compartment is fractal. The fractal nature of this compartment is 


implicit in the operator K, | iF When applied to the corresponding activity concentration 


a; (t) = (1) this operator gives the activity flows to and from compartment F. 


If at least some of the operators K, al | are taken proportional to fractional derivatives of the average 


activities, we obtain a fractional compartment model. The introduction of fractional derivatives does 
not violate mass balance if it is done following a suitable procedure [7]. 

If some of the operators are taken scalar and unsteady, so that some of the average activities are 
multiplied by positive and decreasing time functions to give corresponding flows, we obtain a fractal 
compartment model. In any case, with this general model, the first objective of the present paper is 
attained. 

A simpler version of the proposed three compartments model: 

To simplify the analysis and to be able to construct approximate analytical solutions of tracer kinetics, 
as we already did in the case of the old model, from now on we assume constant compartment 


volumes and substitute the variable volume flows by average constant values. The justification of all 
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these assumptions was given in [4]. Furthermore, we introduce a single operator K [ | and positive 
distribution coefficients @; and > (being Y, + M2 =1) such that: 
(4a) K,,la]= K,,la]=¢, -K[a] (4b) K,,a]=K,,[a]=¢, -K[a] 


Then we obtain simplified equations for tritium kinetics: 


d (2.,+2,) QO, P P, ¢ 
5 A= (Ate vA OP aka te RA 
(Sa) ae V, i() V, 2 V, [ i|+ : [A, ] 
a Gai 4 Sag 5 ply 1 Oo 
(Sb) ria a v, A, V; A, V5 ka, V; Klas] 
a 


(Sc) 


As initial condition we assume an instantaneous and spatially uniform injection of an activity A) of 
tracer in the NE compartment, so A, (0) = Ay In the SW and fractal compartments the initial 


conditions are A (0) =0A (0) = 0 respectively. 

DERIVATION OF A REDUCED ORDER MODEL WITH A’ FRACTAL 
COMPARTMENT AND ASYMPTOTIC ANALYSIS OF TRITIUM KINETICS 

The old model as a reduced order model of the three compartment model: 


If the fractal compartment is neglected we obtain the kinetics already studied in the previous paper 
(old model) [4]. 


Qi. 
6 A A G * 
(6a) i(t)= Ap V5 


jo ,|-Lewpl-Lasl) 


Weal 
(6b) 1l)= | |{H G2 exp s}-(Ji|-92 enh zal 
(7) Ai)= a As0= as oF ol l2,]-2] aoe ral 
Here |A,|<|A,| and (8a) pa [+la,|n@a*@a)., 9 (8b) abel 
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As shown previously [4], the following numerical estimations are obtained (using the hour h as time 


unit): volume V; of SW compartment 3,223 x 10’ m*; volume V2 of NE compartment 1,638 x 107 


m?; average outflow from the SW compartment to coastal waters 2,10 x 10° m/*/h; average 
& cl Pp g 


exchange flows Q,,, Q,, between SW and NE compartments, respectively 8,00x10° m*/h and 1,00x10° 
m°7/h; average fresh water inflow OQ 2 to the NE compartment, 2,00x10° m*/h; average inflow 01. 


from coastal waters to the SW compartment, 1,90x10° m*/h; 


2,|=0,03 h"; 


A,|=0,12 A 
So, there are two time scales: a fast scale |A Ale 8h and aslow scale a.) = 33h 


After 3- la.) = 24h the kinetics in the SW and NE compartments is dominated by the slow time 


scale and both activities A(t) and A,(t) fade in unison with the total activity A(t) in the bay, 


proportionally toexp|- l2,|-¢]. If we define (9a) 7, = and (9b) 7, =l—y, the tritium 


lA | V; 
activities of the SW and NE compartments asymptotically verify: 
(10a) A(t) = x, - A(t) (10b) A,(t) x, - A(t) 
A new reduced order model from the three compartment model: 


Substituting equations (10a) and (10b) in equations (5), adding the first two equations and introducing 
a compartment R whose activity Ag = A, +A, is the sum of the activities in the SW and NE 


compartments, and putting 


1 
(1) A,=|4,| 2) Vp-=V, (13) aoe a ee a (14) A, =A, 


a new two compartments reduced order model follows, including a classical compartment R 


(coincident with the whole bay in the old model) and a fractal compartment F: 


d es 
(15a) one Aw An BAR| + 


1 A 
FRA ERA 


d ier oa 
(15b) arty, BlAgl- 7 Al4e 
R F 
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An approximate solution obtained by perturbation theory: 


If k,is a representative kinetic scalar characteristic of the operator K , we introduce the normalized 


fs ee 1 
operator K,| l= i .K [ iF a representative volume V, such that— =—+—,, and define the 
0 R F 


1 


dimensionless parameter desi! Due to the slowness of mass exchange between the 
Ar Vo (% 
ky 


combination of the SW and NE compartments (compartment R) and the fractal compartment F, we can 


eee 1 
expect that its characteristic time scale —* be small relative to — , and € should be small relative to 1. 
0 R 


This smallness opens the possibility of applying perturbation theory in order to obtain approximate 


solutions describing last stages of the process of tritium exit from the bay. To do this, we introduce 


. k 
two new time scales, an inner (fast) one 7, =A2,-t and an outer (slow) onet, = —*-f. Then, taking 
0 


into account that K,| | is a linear operator and following the standard procedures in perturbation 


theory [8], from equations (15) we obtain two set of equations, two sets of initial conditions and two 
sets of approximate solutions, the so called inner and outer solutions in perturbation theory. Adding 


these internal and external solutions, applying a process of renormalization, and returning to the 


original time, it is possible to construct uniform approximations to the activities A A(t) in compartment 


R, combination of the SW and NE compartments, and A -(t) in the fractal compartment: 


10 
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The function Y (t ) that appears in these formulae is the solution of the initial value problem: 


(17) <“H) ae KW (0). w/(0) =1 (fan operator K is applied to a function f (t) we represent 
t 


F 
the resulting function indistinctly as K [{ of (r)}] or as K [ ve \t)). In a framework of linear fractional 


compartments analysis K[ ]=k, -D? [ ] is proportional to a fractional derivative of order P with 
Fa 


0 < p <1 [7]. Then, if the fractional derivative is taken as Riemann-Liouville [9]: 


(18a) pr)= Bl ky 7 (2) (18b) K{{e*"}]= (zy (ee ee) 


Ve T 
Here E, (z) and EF, B (z)are the one parameter and two parameters Mittag-Leffler functions [9]. The 
function E,, (- | is a monotonically decreasing function of ¢. Plots of this function for 
a =0.25,0.50,0.75,and 1 can be found in reference [9]. When a =1, E, (- !)= exp[—t] 
In a framework of linear fractal compartments analysis [10], K ly |(r)= k(t)- w(t) being I [ |the 
identity operator and k (t ) a positive bounded monotonically decreasing function that tends towards 


zero whent > +00: K[ J=k(t)-7[ | Then: 
AS Heyes faa) 


When the fractal compartment is neglected, all the terms in which the operator K appears vanish. 
Then, from (16a) it follows that A, (1) = Ap (0)-e 7" . But in this case A A(t) must be equal to the 
expression of the asymptotic total tritium activity A(t) of the old model, which is (from equation (7)) 


| 


zA,- . —|A,|-t]. As2,=|A 
A(t) Ay Gay \2.,| t| SH% A, 


injected tritium activity A, thus: 


, the initial value A (0) is related with the total 


___ bl 


(20) A ,(0)= ea) Ao 
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With|Z ,|~0,03A7, |A,|=012 A tand Ay = 100Ci we obtain: A ,(0)=133,33Ci 


_ Aglt)+ Ae(t) 


From the definition of the remnant function for the reduced order model R,,juced (t) = x 0) 
R 


and from equations (16), the following formula for the renewal time T) jeduceq Can be derived: 


1 ¢ 1 1 
— =| Rianeea (t) dt —: = + 
(21) i aaa ea (ea 
Ay \% Vy V, 


1 
If the fractal compartment is neglected, from (21) it follows that the renewal time is — , which is the 
R 


renewal time of compartment R (the combined NE and SW compartments) with only an efflux of 


1 
THO towards the coastal waters. The first term in (21) is smaller than —— because during the first 
R 


stage of the kinetics described by the reduced order model, a net tracer flow from R to F adds to the 
efflux towards the coastal waters, so as far as this is sustained, the water renewal is a little faster than 
in the absence of the fractal compartment. However the second term in (21) is due to the effect of the 
flow reversal between compartments R and F, which occurs in the last stage of tritium kinetics, and 
increases its renewal time. 

With the approximate analytical solutions (16 a), (16 b) and (21), the second objective of the present 


paper is attained. 


FRACTAL ANALYTICAL SOLUTION, RENEWAL TIME AND ASYMPTOTIC BEHAVIOR OF 


TRITIUM REMNANT FUNCTION 


Now, let us suppose that the operator K corresponds to a fractal compartments approach, being 


unsteady and multiplicative: K [ ]=k(t)- i [ | In fractal compartments analysis, to always have bounded 


values of the time dependent and decreasing kinetic coefficients it is usually assumed [5]: 


12 
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—P 
(22) H(t)= hy (1+) O<p<l k,>0 
From (22) it follows: 
(23) f° k(w)-e%" «du =(ky:T)-(Ag-T)?" -€%*7 -T(1- ps4q-T) 


+00 
(Here T(a;x)= fe“ -u“"-du (a>0) isan incomplete Gamma function [11].) 


is a stretched exponential, so: 


hte (EY eyo) kurt (4) 
(24) Hope ty (7) os} cnr) | 


Asymptotically whent — +00, taking into account equations (23) and (24), from the approximate 


solutions (16) it follows: 


ita ae 
(25 % A,(t) : iy (k, T) (Ag T) r(I —P> Ap {E}{ }( Jen el (" 
A,(0) 1+ --(ko-7)-(ay-7)"-e*? TUL “Heda 7) a Vp J\t (l-p) \7 


(25b) Arlt) f 1 
Ap(0) (Ve 


(ko -T)-(Ag-T)?" #7" T= pip 1)} Kot ()"] 


(i-p) \T 
From the asymptotic formulae (25 a) we obtain this asymptotic formula for the tail of the logarithm of 


activity in compartment R: 


~ —p.- J = ko f . f es 
(26) log, A g(t) = (const) p og,( 4 (l- p) ) 


T 


Tentative parameter estimations in this case are (the details of the estimation process will be 


considered elsewhere): p 0.3 T=70h — and - =0.01h" 


F 
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Aj) 4 A(t) 
Whent — + we see, from (25 a) and (25 b), that is dominant over . AS consequence 
A,(0) A,(0) 
A ,(t) 


the tritium remnant function asymptotically vanishes as A,(0)’ that is as a stretched exponential: 
R 


a rb)-{(*7)-0, TY" T(1- pA, 7) ex : fe (es) 


R 


With the asymptotic expression for tritium ’s remnant function (27), our third objective is attained. 


DISCUSSION AND CONCLUSIONS 


Figure 4 shows a logarithmic plot of remnant total tritium activity as a function of time. 


— 


Remnant percent of total injected tracer activity 
eel 


Time since injection (hours) 


Figure 4: Logarithmic plot of tritium activity in the bay (per cent values) as a function of time (in hours, 


with zero corresponding to the injection time). 


A smooth curve was adjusted to measured activities. It has three sectors. 

In the first sector, comprised between the time of the injection and approximately 30 hours, the curve 
begins with a horizontal tangent and is concave down on this time interval. The kinetics is dominated 
by the transference of THO from the NE to the SW compartments. 

In the second sector, comprised approximately between 30 and 70 hours, the logarithmic curve is 


almost straight: the activities in NE and SW compartments fade almost exponentially in unison. 


14 
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In the third and last sector of the curve, it is concave up on the time interval comprised approximately 
between 70 and 100 hours. This reflects the dominance of the transference of THO from the fractal 
compartment to the compartment R. Formula (26) predicts this last concavity. 

There are important limitations to parameters estimation from available data, related with 
measurement uncertainty and low number of measurement data. To diminish this uncertainty, more 
and better measurements are needed, mainly after 100 hours. However the tritium measurements were 
done without electrolytic enrichment of the samples. Corrections were applied for the salinity effects 
on the measured samples of THO [1]. After 100 hours, the uncertainty in the measurements of tritium 
activity in some of the 20 cells that were defined in the bay, made it difficult to estimate the remnant 
total activity. So there was a truncation in the data of total tritium activity in the bay, due to the above 
mentioned reason, closely related with the limitations of the methods of measurement. 


With the old model a mean residence time of 37 h was estimated [4]. This was done applying the 


T. ee ree 
A, |A,| 


formula: Equation (21) modifies this prediction by adding a positive contribution 


to the average renewal time due to the net efflux from F to R. However, the other term in (21) tends to 


diminish slightly the contribution of — to the renewal time, due to the previous net influx from R to 
R 


F. Making a balance between these two opposite tendencies, it is expected an increase in the time of 
residence relative to the value predicted by the old model (37 hours). With the available data of tritium 
activities, additional information is needed to achieve reasonable error bounds that allow a suitable 
estimation of parameters. This is necessary to clarify the importance of the correction due to the effect 


of the fractal compartment. 


A robust mathematical model based on compartments theory, intended to discriminate the tritium 
kinetics in each one of the measured 20 cells, such as the model structure sketched in the discussion 
that appears in reference [4], is not feasible in our case, due to the quality of available experimental 
data and its limited number. The more parameters are involved, the more difficult is the statistical 
demonstration that the model does not in fact describe data. This difficulty appears also in the three 


compartments model if additional information is not used. 


15 
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However, it seems that the construction of a mathematical model with a smaller number of meaningful 
compartments, defined by meaningful additional information about bottom topography and sediments, 
stratifications and velocity fields (from measurements and CFD simulations) in different scenarios, 
like the present one with a fractal compartment, could be useful to analyze experimental data of tracer 
kinetics in other, analogous cases of complex hydrological systems, and to compare with the results of 


purely numerical approaches to these systems. 


The derivation of approximate analytical formulae presented here for a two compartments reduced 
order model was done in a framework that embraces both approaches to compartment analysis: fractal 
and fractional. It remains to investigate the details of the kinetics that predict these two types of 


compartmental analysis and compare them with the available experimental results. 


In sum, the introduction of the fractal compartment allows an explanation to the asymptotic behavior 
of tritium’s remnant function, behavior suggested by the experiments done in Montevideo’s Bay. The 
obtained asymptotic formulae allow the description of non-exponential wash out kinetics from the bay, 
related with the trap effect characteristic of the fractal compartment. This could be of interest to follow 
and assess the impact of the liberation in coastal waters of toxic chemical substances or dangerous 


microbial populations, after a sudden contamination accident in the harbor. 
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